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ABSTRACT 

Using data drawn from the DEEP2 and DEEP3 Galaxy Rcdshift Surveys, we investigate the re- 
lationship between the environment and the structure of galaxies residing on the red sequence at 
intermediate redshift. Within the massive (10 < log]^Q(M^,//i~^ Mq) < 11) early-type population at 
0.4 < z < 1.2, we find a significant correlation between local galaxy overdensity (or environment) 
and galaxy size, such that early-type systems in higher-density regions tend to have larger effective 
radii (by ~ 0.5 kpc or 25% larger) than their counterparts of equal stellar mass and Sersic in- 
dex in lower-density environments. This observed size-density relation is consistent with a model of 
galaxy formation in which the evolution of early-type systems at z < 2 is accelerated in high-density 
environments such as groups and clusters and in which dry, minor mergers (versus mechanisms such 
as quasar feedback) play a central role in the structural evolution of the massive, early-type galaxy 
population. 

Subject headings: galaxicsistatistics, galaxies:fundamental parameters, galaxies :high-rcdshift, galax- 
icsiformation, galaxies:evolution, large-scale structure of universe 



1. INTRODUCTION 

Recent observations of galaxies at intermediate 
redshift (z ^ 2) have identified a significant popu- 
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lation of massive (~ 10 M0), quiescent, early-type 
systems with old, metal-rich stellar populations and 
remarkably sm all sizes relative t o their l ocal coun- 
terparts (e.g., IDaddi et all 120051 : [Labbc ct HI 120051: 
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Tavlor et all |2010[) . The stellar densities of these 
massive, intermediate-redshift galaxies (as measured 
within one effective radius, rg) are typically two orders 
of magnitude greater than quiescent ellipticals of the 
same mass at z ~ 0.1. Within the central 1 kpc 
(physical), however, the densities of early-types at z 2 
are found to only exceed l ocal m e asurements by a factor 
of 2-3 (IBezanson et all 120091: IHopkins et all I2009al: 
Ivan Dokkum et al.l 120101) . Altogether, the observations 
suggest that there is significant evolution in the size of 
massive ellipticals over the past 10 Gyr, likely proceeding 
in an inside-out manner, without the addition of much 
stellar mass. 

Several physical processes have been proposed to ex- 
plain this strong size evolution within the massive, 
early-type population at z < 2. In particular, gas- 
poor, collisionless (or "dry") minor mergers are often 
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). However, a variety of alternative 



mechanisms have also been proposed, including scenarios 
in which the observed structural evolution may be driven 
by secular processes such as adiabatic expansion result- 
ing from stellar mass los s and / or strong AG N- fueled 
feedback (iFaii et al.l [200l l20Tot IDamianov et al. 2003 
IHopkins et al.l l2010bl lal: see also iNipoti et al rf2009l and 
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I Williams et aI1[20lol) . 

One way to possibly discriminate between these sce- 
narios (minor mergers versus secular processes) is by 
quantifying the role of environment in the structural 
evolution of the massive galaxy population. While 
secular processes are largely independent of environ- 
ment and quasars are not p referentially foun d in over- 
dense regions at 2; 1 (jCoil et alj l2007t) . mergers 
are more common in higher-density environments such 
as galaxy groups fCavalicrc ct al.' '1992"; 'Mcintosh ct ajj 
200&: .Wetzel et al. 20 08: Fakhouri & Ma 2009; Lin ct ^ 
2010HDarg et al.l 120101) . Thus, if mergers are the domi- 
nant mechanism by which the sizes of massive early-types 
evolve at z < 2, we should expect to find a variation in 
the structural properties of galaxies as a function of en- 
vironment at z ~ 1. 

To test this, we use data drawn from the DEEP2 
and DEEP3 Galaxy Redshift Surveys (iDavis et"all l200l 
iNewman et al.l 120121 : iCooper et al.l l20lH) to investigate 
the correlation between the local overdensity of galax- 
ies (which we generally refer to as "environment") and 
the sizes of massive galaxies on the red sequence at in- 
termediate redshift. In Section [21 we describe our data 
set, with results and discussion presented in Sections |3] 
and [U respectively. Throughout, we employ a ACDM 
cosmology with w = — 1, = 0.3, Q,a = 0.7, and a 
Hubble parameter of Hq = 100 h km s^^ Mpc^^, unless 
otherwise noted. A ll magnitudes are on the AB system 
(|Oke fc Gunnl[T98l . 

2. DATA 

To characterize both the environment and the struc- 
ture of galaxies accurately requires spectroscopic obser- 
vations (or deep, m ulti-band photometric observations, 
ICooper et al.l I200"5[) as well as high-resolution imaging 
across a sizable area of sky. Given the limitations of 
ground-based adaptive-optics observations, the latter is 
only possible at intermediate redshift via space-based ob- 
servations (e.g., with HST). Among the fields covered by 
deep, multi-band imaging with HST, the Extended Groth 
Strip (EGS) is by far the most complete with regard 
to spectroscopic coverage at intermediate redshift. The 
EGS is one of fou r fields surveyed by the pEEP2 Galaxy 
Reds hift Survey (jPavis et al.ll2003ll2007HNewman et"al] 
I2012D . yielding high-precision (cr^ ~ 30 km s~^) se- 
cure redshifts for 11,701 sources at 0.2 < z < 1.4 over 
roughly 0.5 square degrees in the EGS. Building upon 
the DEEP2 spectroscopic sample, th e recently-completed 
DEEP3 Galaxy Redshift Survey (jCooper et al.l l20TTt 
Cooper et al., in prep) has brought the target sampling 
rate to ~90% at Rab < 24.1 over the central 0.25 square 
degrees of the EGS — th e portion of the field imaged 
by H ST/ ACS (see Fig. [Tl IDavis et all [20071: ILoTz et al.l 
[2OO8I) . 

Among the current generation of deep spectroscopic 
redshift surveys at z ^ 1, the combination of the DEEP2 
and DEEP3 spectroscopic datascts provides the largest 
sample of accurate spectroscopic redshifts, the highest- 
precision vel£cit;j;Jriforni^^ and the highest sampli ng 
density ([Newman et al.|[2"012l : Cooper et al., in prep)["1 

^'^ Note that the samphng density for a survey is defined to be 
the number of galaxies with an accurate redshift measurement per 
unit of comoving volume and not the number of galaxies targeted 



Combined with the relatively wide area imaged with 
HST/ ACS (an area > 2x large r than that surveyed a s 
part of the GOODS program, iGiavalisco et al.l |20M) . 
these attributes make the EGS one of the best-suited 
fields in which to study the relationship between envi- 
ronment and galaxy structure at z '--^ 1. In this paper, 
we utilize a parent sample of 11, 493 galaxies drawn from 
the joint DEEP2/DEEP3 dataset in the EGS with secure 
redshifts (quality Q = 3 o r 4 as defined bv IDavis et al.l 
120071: INewman et al.|[20Tl) in the range 0.4 < z < 1.2. 

2.1. Rest-frame Colors, Luminosities, and Stellar 
Masses 

For each galaxy in the DEEP2/DEEP3 sample, rest- 
frame U — B colors and absolute B-haiid magni- 
tudes, Mb, are calculated from CFHT BRI photom- 
etry ([Coil et all [200I using the JC-co rrection proce- 
dure described in iWillmer et al.1 ([2006D . For a sub- 
set of the galaxy catalog, stellar masses are calcu- 
lated by fitting spectral energy distributions (SEDs) to 
WIRC/Palomar J- and Kg-hand photometry in conjunc- 
tion with the DEEP2 BRI data, according to t he pre- 
scriptions described bv lBundv et all ([20051I2006D . How- 
ever, the near-infrared photometry, collected as part of 
the Palomar Observa tory Wide-field Infrared (POWIR, 
iConsehce et"an [2001) survey, does not cover the entire 
DEEP2/DEEP3 survey area, and often faint blue galax- 
ies at the high-z end of the DEEP2 redshift range are 
not detected in Kg ■ Because of th ese two effects, the 
stellar masses of IBundv et al.l ([20061 ) have been used to 
calibrate stellar mass estimates for the full DEEP2 sam- 
ple that are based on rest-frame Mb and B — V values 
derived from t he DEEP2 data in conjunction with the 
expressions of [Bell et all ([2003D . which relate mass-to- 
light ratio to optical color. W e empirically c orrect these 
stellar mass estimates to the IBundv et al.l ([2006() mea- 
surements by accounting f or a mild color and redshift de- 
pendence ([Lin et al.ll2007[) : where they overlap, the two 
stellar masses have an rms difference of approximately 
0.3 dex after this calibration. 

2.2. Local Galaxy Overdensity 

To characterize the local environment, we compute 
the projected third-nearest-neighbor surface density (S3) 
about each galaxy in the joint DEEP2/DEEP3 sam- 
ple, where the surface density depends on the projected 
distance to the third- nearest neighbor, i?p,3, as S3 = 
3/(7r£>p 3). In computing S3, a velocity window of ±1250 

km s~^ is utilized to exclude foreground and background 
galaxies along the line of sight. Varying the width of 
this velocity window (e.g., using ±1000-2000 km s^^) 
or tracing environment according to the projected dis- 
tance to the fifth-nearest neighbor has no signifi c ant ef - 
fect on our results. In the tests of ICooper et al.l ([2005D . 
this projected n*^-nearest-neighbor environment estima- 
tor proved to be the most robust indicator of local galaxy 
density for the DEEP 2 survey. 

To correct for the redshift dependence of the DEEP2 
and DEEP3 sampling rates, each surface density value 
is divided by the median S3 of galaxies at that redshift 
within a window of Az = 0.04; correcting the measured 

down to an arbitrary magnitude hmit. 
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surface densities in this manner converts the S3 values 
into measures of overdensity relative to the median den- 
sity (given by the notation 1+83 here) and effectively 
accounts for the red shift variations in the selection rate 
(|Cooper et alJ 120051 ). Finally, to minimize the effects of 
edges and holes in the survey geometry, we exclude all 
galaxies within 1 h"^ comoving Mpc of the DEEPS sur- 
vey boundary (see Fig. [1]), reducing our sample to 7, 257 
galaxies in the redshift range 0.4 < z < 1.2. 
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Fig. 1. — DEEP2 and DEEPS spectroscopic coverage in the Ex- 
tended Groth Strip. The black and green points denote the spec- 
troscopic targets of the two surveys, with the area covered by 
HST/ACS imaging highUghtcd in orange. The dashed magenta 
line denotes the edge of the DEEPS survey and thus the area over 
which environments were computed for this work. The areas to be 
imaged with HST/ACS and W FCS-IR as part of the CANDELS 
Multi-Cycle Treasury Program I IGrogin et al.|[2blll) are delineated 
by the cyan and red outlines, respectively. Finally, the location of 
the fo ur CFHT pointings (numbered 11-14) imaged bv lCoil et al.l 
l|2004l ) are denoted by the black dashed lines. 



2.3. Sersic Indices and Sizes 

To quantify the sizes of the DEEP2 and DEEPS galax- 
ies, we utilize morphological measurements extracted 
from the Advanced Camera for Surveys General Cata- 
log (ACS-GC, Griffith et al. 2012, in prep). The ACS- 
GC analyzed the HST/ACS Vf60 6W and /f814W imag ing 
in the EGS using GALAPAGOS ()Haufiler et al.llMll) . a 
semi-automated tool for measuring sizes and spatial pro- 
files via the parametric fitting code GALFIT (jPeng et al.l 
l2002l[2010f ). To determine the galaxy profile shape, each 
radial profile was fit using a simple Sersic measurement 
(|Sersidll968D of the fo rm 



T,{r) = Ee exp[-K((r/re)~ 



1)] 



(1) 



where is the effective radius of the galaxy, Se is the 
surface brightness at re, n is the power-law index, and k 
is coupled to n such that half of the total flux is always 
within Tg. 

Here, we employ the profile fits to the HST/ACS 
^F8i4W imaging for all sources, independent of galaxy 
redshift. At 2 < 1.2, the /f814W passband samples 
the rest-frame optical (A > S700A), which minizes mor- 
phological biases associated with observations made in 
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Fig. 2. — Top: the rest-frame U — B versus Mb color-magnitude 
distribution for the DEEP2 and DEEPS galaxies in the spectro- 
scopic sample within the redshift range 0.4 < z < 1.2. The 
four solid black lines show the effective absolute-magnitude limit 
of the survey at z = 0.7, 0.9, 1.05, and 1.2, while the solid and 
dashed red lines show lines of constant stellar mass corresponding 
to logiQ(M*//i~2 M0) = 10 and 10.5, respectively. The dashed 
blue/red horizontal line shows the division betwe en the red se- 
quenc e and the blue cloud as given by Equation 19 of lWillmer et all 
II2OO6I ). Bottom: the combined distribution of the 16,676 unique, 
secure (Q = -1, S, 4) redshifts measured by DEEP2 and DEEPS 
within the Extended Groth Strip field (black histogram). The 
red and blue histograms show the corresponding distributions for 
the DEEP2 and DEEPS surveys independently. Recall that the 
DEEP2 survey covers roughly a factor of 2 more area than the 
DEEPS survey (see Fig.[T]|. 



the rest-frame ultraviolet where galaxies typically ex- 
hibit more irregular morphologies. The impact (or lack 
thereof) of any morphological if-correction is addressed 
further in iJ21 Throughout our analysis, all sizes (rg) 
have been converted to physical kpc, according to the 
DEEP2/DEEPS spectroscopic redshift and assuming a 
Hubble parameter oi h = 1. Finally, note that the mul- 
tidrizzlcd HST/ACS images, from which structural prop- 
erties were measured, have a pixel scale of 0.03" per pixel 
and a point-spread function (PSF) of 0.12" FWHM; from 
z = 0.4 to z = 1.2. the spatial resolution therefore varies 
from -0.5 kpc per PSF FWHM to -0.95 h'^ kpc 
per PSF FWHM. 

2.4. Sample Selection 

To investigate the relationship between galaxy struc- 
ture and environment amongst the high-mass por- 
tion of the red sequence, we define a subsample of 
DEEP2/DEEPS galaxies at 0.4 < z < 1.2 with stellar 
mass in the range 10 < log]^o(^*/^~^ ^0) < the 
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red sequence (i.e., rest- frame color oiU — B > 1)13 and 
with robust environment and morphology measurements 
(i.e., away from a survey edge and with cr„ < 0.75)13 
The median redshift for this subsample of 623 galaxies is 
0.76 and the median stellar mass is log]^Q(M^//i^^ M©) ~ 
10.6. 

In Figure [U we show the redshift distribution for all 
sources in the EGS with a secure redshift {Q ~ —1, 3, 4) 
in the joint DEEP2/DEEP3 sample alongside the color- 
magnitude distribution for all galaxies at 0.4 < z < 1.2, 
with lines of constant stellar mass overlaid and with lines 
illustrating the survey magnitude limit at several discrete 
redshift values. Note that we restrict our primary sub- 
sample to a redshift range over which the DEEP 2 and 
DEEP3 selection function is relatively flat. However, 
over this somewhat broad redshift range the sample is 
incomplete at the adopted mass limit. For example, at 
z = 0.9 the Rab = 24.1 magnitude limit of DEEP2 in- 
cludes all galaxies with stellar mass > 10^"'^ M^,//i~^ Mq 
independent of color, but preferentially misses red galax- 
ies at lower masses (see Fig. [2]). This incompleteness in 
the galaxy population is addressed in more detail in 

3. ANALYSIS 

In order to study the relationship between galaxy prop- 
erties and environment at fixed stellar mass, as we aim 
to do here, the galaxy sample under study is often re- 
stricted to a narrow range in stellar mass such that cor- 
relations between stellar mass and environment are negli- 
gible. However, at intermediate redshift, sample sizes are 
generally limited in number such that using a particularly 
narrow stellar mass range (e.g., 0.1-0.2 dex in width) 
significantly reduces the statistical power of the sample. 
For this reason, broader stellar mass bins (e.g., ~ 0.5 
dex) are commonly employed. However, if the shape of 
the stellar mass function depends on environment (as 
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the typical stellar mass within a broad mass bin may 
differ significantly from one density regime to another. 
Such an effect would clearly impact the ability to study 
the relationship between galaxy size and environment at 
fixed stellar mass. 

For this reason, we instead select those galaxies within 
the top 15% of the overdensity distribution for all red 
galaxies at 10 < logio(M*//i~2 < n and 0.4 < 

z < 1.2 (a subsample of 93 galaxies), and from the corre- 
sponding bottom 50% of the overdensity distribution we 
randomly draw 1000 subsamples (each composed of 93 
galaxies) so as to match the joint redshift, stellar mass, 
and Scrsic index distributions of the galaxies in the high- 
density subsamplcPI The average environment for the 
high-density subset is log]^g(l -I- S3) = 1.31, with an in- 

We adopt this simplified color-cut to be more restrictive 
at the faint end of the red sequence, where dusty star-forming 
galaxies are more prominent IILotz et al.ll2009) . Ho wever, using a 
lumin osity-dependent color-cut (e.g., Equation 19 of lWillmer et all 
1200611 yields no significant changes in our results. 

^® Limiting the sample to those sources with (T„ < 0.75 excludes 
very few (only 2 out of 625) objects. Removing this selection cri- 
terion or making it more restrictive (e.g., cr„ < 0.5) yields no sig- 
nificant changes in our results. 

"^^ Selecting the top 10% or 20% of the environment distribution 
yields similar results. 



terquartile (25%-75%) range of 1.15-1.41, while the low- 
density subsample is biased to considerably less-dense 
environs with an average overdensity of log]^Q(l -I- Js) = 
—0.12 an d an interquartile range of —0.30-0.12. As 
shown bv lCooper et al.l ()2006[ ). the high-density subsam- 
ple is comprised largely of group members, while the low- 
den sity population is dominated by "field " galaxies (see 
also lCooper eraII[2007l:[Gerke et al.l[2007l) . 

By matching in redshift, we remove the projection of 
any possible residual correlation between our environ- 
ment measurements and redshift (in concert with the 
known redshift dependence of the survey's stellar-mass 
limit) onto the observed size-environment relation. In 
addition, matching according to redshift alleviates any 
possible impact from morphological if-corrections as- 
sociated with measuring all structural parameters in 
the HST/ ACS /f814w passband. Finally, recognizing 
the correlations between environment and parameters 
such as color, star-formation rate, and mo rphology (i.e., 
early - versus late-type) at z ^ 1 (e.g.. iCooper et al.l 
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al.ll2007h . we also force the Scrsic indices 
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of the low-density subsample to match those of the high- 
density population, which controls for the contribution 
of dusty disk galaxies to our red-sequence population. 
Matching based on rest-frame color, in lieu of Sersic in- 
dex, would confuse reddened disk galaxies with red early- 
type systems. 

Members of the low-density subsample are drawn ran- 
domly from within a three-dimensional window with di- 
mensions of |Az| < 0.04, |Alogio(M^//i-2 Mq)\ < 0.2, 
and I An| < 1.5 of a randomly-selected object in the high- 
density subsample. Varying the size of this window by 
factors of a few in each dimension has no significant effect 
on our results. Given the random nature of the matching, 
some objects are repeated in the low-density subsamples. 
However, for each subsample of 93 galaxies, > 70% of the 
galaxies are unique; requiring all members of a subsam - 
ple to be unique would skew the statistics (|Efronlll981|) . 
By matching our high- and low-density subsamples in 
stellar mass, redshift, as well as Sersic index, we are able 
to effectively study the correlation between galaxy size 
and environment at fixed stellar mass. 

To test whether our high-density subsample and 
the random low-density subsamples arc consistent with 
being drawn from the same underlying stellar mass 
distribution, we apply two non-parametric (i.e., in- 
dependent of Gaussian assumptions) tests, the two- 
sided Kolmogorov-Sm irnov (KS) test (IPress et all 119861: 
IWaU fc JenkiS [20031) and the one-sided Wilcoxon- 
Mann - Whitney (WMW) U test (|Mann fc WhitnevI 
Il947| ). The result of each test is a P- value: the prob- 
ability that a value of the KS or U statistic equal to the 
observed value or more extreme would be obtained, if 
the "null" hypothesis holds that the samples are drawn 
from the same parent distribution. The WMW U test is 
computed by ranking all elements of the two data sets 
together and then comparing the mean (or total) of the 
ranks from each data set. Because it relies on ranks 
rather than observed values, it is highly robust to non- 
Gaussianity. The WMW U test is particularly useful for 
small data sets (e.g., compared to other related tests suc h 
as the chi-square two-sample test. I Wall fc Jenkin j|2003() . 
as we have when selecting galaxies from a narrow stel- 
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lar mass range and in extreme environments, due to its 
insensitivity to outlying data points, its avoidance of bin- 
ning, and its high efHciency. Note that since this test is 
one-sided, possible Pjj values range from to 0.5 (versus 
Pks which ranges from to 1); for a Pjj value below 0.025 
(corresponding closely to 2a for a Gaussian), we can re- 
ject the null hypothesis (that the two samples have the 
same distribution) at greater than 95% significance £3 

In Figure [3l we plot the cumulative distribution of 
stellar masses for the 93 sources in the high-density 
subsample alongside that for the 1000 random subsam- 
ples (each consisting of 93 galaxies) matched in redshift 
but residing in low-density environments. Performing 
a one-sided WMW U test (and a two-sided KS test) 
on the size (re) measurements for the low- and high- 
density populations, we find that the size distribution 
for the galaxies in high-density environments is skewed 
to larger sizes, with a probability of Pjj < 0.01 (and 
^ks < 0.02). Meanwhile, the cumulative stellar mass, 
redshift, and Sersic index distributions for the low- and 
high-density subsamples, shown in the inset of Figure 
[21 arc well-matched with the WMW U test yielding a 
Pu > 0.4. This confirms that our sample-construction 
procedure has yielded sets of galaxies in low- and high- 
density environments whose redshift, mass, and Sersic 
index distributions match closely. While not directly 
matched, the rest-frame color distributions for the two 
samples are also indistinguishable — not a surprising re- 
sult given that the color-density relation shows no sig- 
nificant v ariation across the r e d sequence at a giv en lu- 
minosity (jBlantoEiEaD IIMI ICooEerira^ See 
Table [1] for a complete summary of the probability values 
given by both the WMW U and KS tests. 

The results of the WMW U test are confirmed by a 
comparison of the Hodges-Lehmann (H-L) estimator of 
the mean sizes for the low- and high-density subsamples, 
which differ by 0.54 ± 0.22 kpc. This reinforces the 
conclusion that there is a non-negligible size-environment 
relation on the red sequence at z ^ 0.75. The Hodges- 
Lehmann (H-L) estimator of the mean is given by the 
median value of the mean computed over all p airs of 
galaxies in the sample (jHodges fc Lehmaiinlll963[ ). Like 
taking the median of a distribution, the H-L estimator 
of the mean is robust to outliers, but, unlike the median, 
yields results with scatter (in the Gaussian case) compa- 
rable to the arithmetic mean. Thus, by using the H-L 
estimator of the mean, we gain robustness as in the case 
of the median, but unlike the median, our measurement 
errors are increased by only a few percent. 

In Figure HI we show the distribution of the differ- 
ences between the Hodges-Lehmann estimator of the 
mean size, stellar mass, redshift, Sersic index, and color 
for the high-density subsample relative to that for each 
of the 1000 low-density subsamples, where the me- 
dian difference in the H-L estimate of the mean size 
is Are = 0.559 (as illustrated by the dashed verti- 
cal line) versus A logio(M*//i-2 Mq) = -0.002, Az = 
-0.003, An = 0.057, and A{U - B) = -0.001 for 
stellar mass, redshift, Sersic index, and color, respec- 
tively (see Table [2]). Within the stellar mass range of 

21 Note that the two-sided probability for the WMW U test is 
given by doubling the one-sided probability. Here, we report only 
the one-sided probability. 



10 < logiQ(M^//i~^ Mq) < 11, we find significant evi- 
dence for a correlation between galaxy size and environ- 
ment at z 0.75, such that higher-density regions play 
host to larger galaxies at a given stellar mass on the red 
sequence. 
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Fig. 3. — the cumulative size (re), stellar mass, redshift, Sersic 
index (n), and rest-frame color (U — B) distributions for the 93 
DEEP2 galaxies comprising the top 15% of the environment dis- 
tribution within the stellar mass and redshift ranges of 10 < 
logiQ(M*//i-2 Mq) < 11 and 0.4 < ^ < 1.2 in comparison to the 
corresponding cumulative distributions for the 1000 random galaxy 
subsamples drawn from the lowest 50% of the same environment 
distribution. As discussed in the text, the low-density subsamples, 
which are each composed of 93 galaxies, are selected to have the 
same stellar mass, redshift, and Sersic index distribution as the 
high-density population. However, the size distribution is found to 
be significantly different in the different environment regimes. 

To test the robustness of our results to the particular- 
ities of the sample selection, we repeat the analysis de- 
scribed above for several samples spanning varying red- 
shift and stellar mass regimes. For example, restrict- 
ing the redshift range over which we select galaxies to 
0.7 < z < 1.2, thereby decreasing the size of the sam- 
ple, we again find a statistically significant relationship 
between galaxy structure and environment within our 
adopted stellar mass bin of 10 < log]^Q(Mj,//i^^ ^q) < 
11. For the 64 galaxies in the high-density regime 
(again the highest 15% of the overdensity distribution) 
at 0.7 < z < 1.2, the cumulative distribution of galaxy 
size (re) is skewed towards larger effective radii relative 
to the comparison set of galaxies in low-density environ- 
ments, yielding Pjj < 0.01 and a median difference in the 
H-L estimator of the mean size of Are ~ 0.63 kpc. 

An obvious concern when studying the structure of 
galaxies on the red sequence is the relative contribution of 
early-type systems and dusty disk galaxies to the sample. 
Especially at fainter magnitudes/lower stellar mass, red- 
dened star-forming galaxies comprise a significant por- 
tion of the red galaxy population (e. g., iLotz et 311120081 : 
IBundv et al.l[20iaiCheng eFalllMTh . and as many stud- 
ies of environment and clustering at intermediate redshift 
have shown the star-forming population tends to reside 
in lower-de nsity regions relativ e to their passive counter- 
parts fe.g. iCooper et aL|[20ql [20071 : iCapak et aL| [2007l: 
iCoil et al.1120081: iKovac et al.ll2010[) Even with our pri- 
mary sample selected to be at the massive end of the red 
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within the high-mass segment of the early-type galaxy 
population the size-stellar mass relation varies with en- 
vironment, such that galaxies in high-density regions are 
larger at a given stellar mass. 



-0.5 0.0 0.5 



Fig. 4. — the distribution of differences between the Hodges- 
Lehmann (H-L) estimator of the mean size, stellar mass, red- 
shift, Sersic index, and rest-frame color for the high-density sub- 
sample of 93 red-sequence galaxies at 0.4 < z < 1.2 and 10 < 
logj^o(M*/fc^^ M0) < 11 relative to the corresponding H-L estima- 
tor of the mean for each of the 1000 low-density subsamples. The 
dotted and dashed vertical lines show the median and mean, re- 
spectively, of the distribution of differences between the H-L means 
for each galaxy property; error bars denote the uncertainty in the 
mean as given in Table [2] We find a significant offset in size of 
Are ~ 0.6 kpc (physical), while the difference in mean stel- 
lar mass, redshift, Sersic index, and color of the two samples is 
consistent with zero (by construction for all but color). 

sequence, disk galaxies {n < 2.5) still account for '^25% 
of the population (see Figure It is unlikely, however, 
that a difference in the relative contribution of dusty disk 
galaxies to the high- and low-density populations is driv- 
ing the observed size-density relation since galaxy size 
correlates with Sersic index such that red disk galax- 
ies (systems with n < 2.5) tend to have slightly larger 
(not smaller) measured sizes than red galaxies of slightly 
larger Sersic index. Still, to minimize any potential im- 
pact from late-type systems, we define two subsamples: 
one selected according to n > 2.5 and a second further 
constrained to systems with n > 2.5 and with axis ratios, 
(6/a), greater than 0.4 (see Table [T]). For these subsam- 
ples of galaxies with early-type morphology, we find that 
the correlation between structure and environment per- 
sists, reinforcing the conclusion that early-type systems 
(and not dusty disk galaxies) arc responsible for the ob- 
served size-density relation at fixed stellar mass. 

To illustrate the correlation between size and environ- 
ment in a more physically intuitive manner, in Figure 
O we show the size-stellar mass relations for early-type 
galaxies in overdense and underdense regions. In the 
top portion of Fig. [51 the high- and low-density sam- 
ples are simply defined to be the extreme quartiles of the 
environment distribution for all galaxies with n > 2.5, 
0.7 < z < 1.2, and U ~ B > I. In the middle and 
bottom panel, however, a more controlled comparison is 
made, with the high-density sample selected as the top 
15% of the environment distribution for all galaxies with 
n > 2.5, 0.4 < z < 1.2, and U - B > 1 and the low- 
density sample is comprised of 400 galaxies randomly 
drawn from the bottom 50% of the environment distribu- 
tion to match the z distribution (middle) and the joint n 
and z distribution (bottom) of the high-density sample. 
In agreement with our previous analysis, we find that 
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Fig. 5. — the mean (filled diamond connected by dashed lines) 
and median (open hexagons connected by dotted lines) relation- 
ships between galaxy size (re) and stellar mass in the high-density 
(red lines and symbols) and low-density (blue lines and symbols) 
regimes. The means and medians are computed in bins of stellar 
mass with width A log]^Q(Mj,//i~^ M0) = 0.3. In the top panel, 
the high- and low-density samples are selected as the respective 
extreme 25% of the environment distribution for all galaxies with 
U — B > 1, 0.7 < z < 1.2, and n > 2.5 — with no matching accord- 
ing to n or z. In the middle and bottom panels, the high-density 
sample comprises the top 15% of the environment distribution for 
all galaxies with U — B > I, 0.4 < z < 1.2, and n > 2.5, with the 
low-density sample selected from the bottom 50% of the environ- 
ment distribution so as to have the same z (and n) distributions 
(but not stellar mass). In both instances, we find that the size- 
stellar mass relation for early-type galaxies is systematically offset 
to larger sizes in overdense regions. 

Finally, at low surface brightness levels (i.e., low signal- 
to- noise per pixel in the HST/ ACS imaging), GAL- 
FIT tends to underestimate both galaxy size (rg) and 
Sersic index. While such effects are minimal for ob- 
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jects brighter than the sky background (|Haufiler et alj 
we define a subsample limited to those sources 
with surface brightness in the HST/ ACS /f814W pass- 
band of ^ < 23.5 magnitudes per arcsec^ (A^sky ^ 27.4), 
which excludes 58 of 623 red galaxies at 0.4 < z < 1.2 
and 10 < logio(M*//i^2 Mq) < 11. For our samples 
with closely matched stellar mass, z, and color (i.e., ef- 
fectively matched apparent magnitude), such a cut on 
surface brightness imposes an upper limit on physical 
size, which therefore may impact the observed strength 
of the size-density relation. However, in spite of this 
conservative surface brightness limit, we still find a sig- 
nificant relationship between size and local overdensity 
on the red sequence, with a median difference in the H-L 
estimator of the mean size of Ar^ ^ 0.47 kpc. In 
Table [Hand Tabic [21 we list the results from similar anal- 
yses of several different galaxy samples. When varying 
the Sersic index, stellar mass, surface brightness, axis ra- 
tio, and/or redshift regimes probed, we continue to find 
a significant size-density relation at fixed stellar mass on 
the red sequence. 

4. DISCUSSION 

Previous efforts to study the environmental depen- 
dence of the size-stellar mass relation at z < 2 within 
the early-type galaxy population have been relatively 
few in number and have often found no significant 
trends with local galaxy density. Comparing the mor- 
phologies of massive galaxies in the low-redshift (z = 
0.165) Abell 901/902 supercluster to those of comparable 
field samples selected from the Space Telescope A901/ 2 
Galaxy Evolu tion Survey f STAGES. iGrav et all [2009I) . 
iMaltbv et all (I2010D detect no significant relationship 
between galaxy structure and environment within the 
early-type population. Focusing on galaxy g roups iden- 
tified in the Sloan Digital Sky Survey fSDSS. lYork et al.l 
|2000[) . a less extreme subdivi sion of the environment 
distribution, iGuo et aTl ()2009[ ) also find no significant 
evidence for a correlation between local environment 
and the size or Sersi c index within th e loca l early- type 
popu lation (see also iWeinmann et al.l l2009t iNair et aD 
120101 ). In contrast, previous studies of brightest clus- 
ter galaxies (BCGs) in the local Universe find that 
BCGs tend to be larger than early-types of comparable 
stellar mass in less-dense environs (e.g.. [Bernardi et al.l 



2007t Ivon der Linden et al.lj 2007: 'Dcsrochcs et al."20q 
Liu et all I2008L see also iBovlan-Kolchin et al., ,200 ' 



However, this environmental dependence apparent in the 
BCG population is likely the result of cluster-specific 
mechanisms that drive the formation of this rare sub- 
set of the massive early-type galaxy population and may 
not be indicative of the massive early-type population 
as a whole. We note that our DEEP2/DEEP3 sample 
includes very few (if any) systems that will evolve into 
BCGs at z - 0. 

Beyond the local Universe, iValentinuzzi et all ()2010[ ) 
find no significant variation in the size-stellar mass re- 
lation for massive early- types in clusters at z ~ 0.7 rel- 
ative to that for comparable systems in the field, us- 
ing data drawn from the ES Q Distant Clusters Survey 
(EDisCS. I White et al.l [200l. At yet higher redshift, 
iRettura et al.l ( 2010[) compare the sizes of massive ellip- 
ticals in an X-ray- luminous cluster at z = 1.237 to those 
of correspondingly-massive systems in the field, finding 



no significant variation in size wit h environment . Uti - 
lizing the same data set, however, iCimatti et al.l (|2008f ) 
propose (though without quantifying) a possible corre- 
lation between size and environment similar in nature 
to that found within our DEEP2/DEEP3 sample (that 
is, higher-density regions favoring less-compact galax- 
ies). Nevertheless, using largely the same galaxy sam- 
ples as these two pr evious studies, recent work from 
iRaichoor et~al] (|201lD argues for the opposite trend such 
that galaxies in high-density regions are smaller than 
their field counterparts. The significance of the measured 
correlation between size a nd environment , howe ver, is 
dramatically overstated by IRaichoor et al.l (|2011[ ) , with 
their results actually consistent with no environment de- 
pendencel3 The lack of a significant environment depen- 
dence reported in these previous studies i s likely due to 
i] the sma ller sample sizes emp loyed (the IRettura et aTl 
20101 ) and ICimatti et al.l (|2008| ) analyses included a to- 
tal of 45 ellipticals across all environments), [ii] the use of 
less-precise environment measures (e.g., relying on pho- 
tometric redshifts such that "field" or low-density sam- 
ples can be strongly contaminated by group members), 
and/or [in] differences in the redshift range probed. 

Our results, which show a significant size-density rela- 
tion at fixed stellar mass within the massive, red galaxy 
population, suggest that the structural evolution of mas- 
sive early-type systems occurs preferentially in overdense 
environments (i.e., g roups, given the la c k of m assive clus- 
ters in our sample, iGerke et al]l2005L l2012| ). This en- 
vironmental dependence is in general agreement with a 
model of galaxy formation in which minor, "dry" merg- 
ers are a critical mechanism in driving structural evo- 
lution at z < 2. Moreover, our DEEP2/DEEP3 re- 
sults also suggest that early-types in higher-density re- 
gions evolved structurally prior to their counterparts in 
low-density regions, growing from highly-compact sys- 
tems at z ~ 2-3 to more extended systems at z ~ 0. 
The earlier onset of this evolution in overdense envi- 
ronments is in agree ment with studies of stellar pop- 
ulations locally (e.g., iCooper et al.l l2010al ) as well as 
studies of the color-density relation at intermediate red - 
shift (e.g., ICooper et all [20061 \2(M . IGerke et al.l[2007[) . 
which find that galaxies in high-density environments 
have typically ceased their star formation earlier than 
those in less-d ense environs. Studies of t he Fundamen- 
tal P lane (FP, iDjorgovski fc DavisI [l987l : [Dressier et "aD 
Il987| ) support this picture of accelerated evolution in 
high-density environments, with several analyses find- 
ing that galaxies in high-density regions tend to reach 
the F P more quickly than th o se in low-density re- 
gions (!van Dokkum et al.' '2001': 'Ge bhardt efall 120031: 
iTreu et al. 2005a; Moran ct al. 2005). 

A correlation between size and environment at fixed 
stellar mass within the massive early-type population is 
arguably in conflict with scenarios in which the observed 
size evolution at intermediate redshift is driven by quasar 
feedback, as quasars at z 1 are generally not found 
to reside in overdense environments, especially in rela- 
tion to the early-type galaxy population. Using cross- 
correlation techniques and measurements of local envi- 

22 The average s izes for cluster, group , and field early-types as 
given m Table 1 of lRaichoor et al.l pOTll ^ are all consistent at the 
Icr level. 
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TABLE 1 

Summary of KS and WMW U Probabilities 



Sample 


-^high — dcnsty 


logio(M*) 
Pks Pu 


Pks 


Pu 




n 

Pu 


U ■ 
Pks 


- B 
Pu 


Pks 


Pu 




U.4 < 2 < i.Z 

10 < logio(M*) < 11 


93 


0.718 


0.480 


0.605 


0.428 


0.220 


0.413 


0.542 


0.304 


0.011 


0.006 


0.7 < 2 < 1.2 

10 < logio(M*) < 11 


64 


0.887 


0.421 


0.578 


0.356 


0.562 


0.461 


0.321 


0.280 


0.036 


0.002 


0.4 < 2 < 1.05 
10 < logio(M*) < 11 


86 


0.807 


0.493 


0.740 


0.445 


0.330 


0.395 


0.614 


0.387 


0.042 


0.020 


0.4 < 2 < 1.2 
10.5 < logi(,(M*) < 11 


58 


0.795 


0.494 


0.868 


0.452 


0.220 


0.391 


0.379 


0.246 


0.015 


0.013 


0.7 < 2 < 1.2 
























10 < logio(M*) < 11 


42 


0.729 


0.475 


0.725 


0.313 


0.230 


0.384 


0.301 


0.298 


0.007 


0.002 


n > 2.5 
























0.7 < 2 < 1.2 
























10 < logn,(M*) < 11 


36 


0.711 


0.476 


0.654 


0.377 


0.252 


0.327 


0.649 


0.393 


0.002 


0.001 


n > 2.5, b/a > 0.4 
























0.4 < 2 < 1.2 
























10 < logio(M*) < 11 


85 


0.810 


0.499 


0.525 


0.453 


0.594 


0.411 


0.796 


0.364 


0.059 


0.010 


At < 23.5 



























Note. — For several galaxy samples, we list the results of the KS and WMW U tests (Pks ^nd Pjy, respectively) 
from a comparison of the stellar mass, redshift, Scrsic index, color, and size distributions of the respective low- and 
high-density samples. The P-values indicate the probability that differences in the distribution of the stated quantity 
as large as those observed (or larger) would occur by chance if the two samples shared identical distributions. Recall 
that the P[/ values are one-sided probabilities, while the Pks are two-sided. The number of galaxies in the high-density 
sample (picked to be the top 15% of the environment distribution) is given by A^jjigij. density Note that stellar masses 
and sizes are in units of h~-^ Mq and physical kpc, respectively. The samples listed here are matched in stellar 
mass, redshift, and Scrsic index, by construction. 



TABLE 2 

Summary of Differences between Hodges-Lehmann Estimates of the Mean 



Sample 


logio(M*) 


2 


Ahl 

n 


U - B 


re 


median of A 

logl0(M*) 2 


HL distribution 
n U - B 


re 




0.4 < 2 < 1.2 

10 < logio(MO < 11 


-0.003 
±0.030 


-0.004 
±0.013 


0.047 
±0.248 


0.001 
±0.011 


0.543 
±0.220 


-0.002 


-0.003 


0.057 


-0.001 


0.559 


0.7 < 2 < 1.2 

10 < logio(M,) < 11 


-0.022 
±0.026 


-0.004 
±0.008 


0.023 
±0.272 


-0.008 
±0.012 


0.688 
±0.251 


-0.005 


-0.003 


0.075 


-0.010 


0.632 


0.4 < 2 < 1.05 
10 < logio(M*) < 11 


0.010 
±0.024 


-0.005 
±0.011 


0.085 
±0.226 


-0.005 
±0.011 


0.393 
±0.209 


-0.003 


-0.002 


0.078 


0.001 


0.449 


0.4 < 2 < 1.2 
10.5 < logio(M*) < 11 


-0.008 
±0.017 


0.001 
±0.015 


0.150 
±0.285 


0.010 
±0.012 


0.619 
±0.274 


-0.003 


-0.002 


0.085 


0.012 


0.625 


0.7 < 2 < 1.2 
10 < logi(,(M*) < 11 
n > 2.5 


-0.002 
±0.030 


-0.002 
±0.009 


0.145 
±0.277 


-0.010 
±0.013 


0.802 
±0.239 


-0.002 


-0.003 


0.090 


-0.006 


0.826 


0.7 < 2 < 1.2 
10 < logio(M,) < 11 
n > 2.5,b/a > 0.4 


-0.003 
±0.030 


-0.005 
±0.016 


-0.128 
±0.301 


-0.006 
±0.013 


0.909 
±0.273 


-0.001 


-0.006 


0.175 


0.000 


0.981 


0.7 < 2 < 1.2 

10 < logio(M*) < 11 
At < 23.5 


0.000 
±0.024 


0.001 
±0.011 


0.083 
±0.214 


-0.001 
±0.010 


0.426 
±0.224 


0.000 


0.000 


0.050 


0.000 


0.471 





Note. — For the same galaxy samples listed in Table[T] we list the mean and median of the distribution of the differences 
in the Hodges-Lehmann estimator of the mean stellar mass, redshift, Sersic index, color, and size for the high-density samples 
relative to that for each of the low-density samples (see Figure|4]l. Note that stellar masses and sizes are in units of Mq 
and physical kpc, respectively. Finally, note that the samples listed here are matched in stellar mass, redshift, and 
Sersic index, by construction. 



ronment analogous to those presented herein, iCoil et all 
()2007D find that quasars at z ~ 1 cluster like blue galax- 
ies, such that they favor regions of average galaxy den- 
sity (logiQ(l + S3) ~ 0). While these results suggest that 
quasars are unlikely to be responsible for the larger sizes 
of early-types in high-density regio ns, it should be n oted 
that the clustering measurements of lCoil et all (l2007l) are 
roughl y l-2g lower than th at of si milar studies. For ex- 
ample. ICroom et ahl ()2005l see also iPorciani et ai1l2004l : 



IGrazian et "alll2004l: IMvers et al']|2006D find that quasars 
at 0.3 < z < 2 . 2 hav e a bias > 2a higher than th at found 
bv ICoil et al.1 (l2007l) . while ISerber eFall (l2006l) find an 
excess of ~ L* galaxies on 25 kpc to 1 Mpc {h = 0.7) 
projected comoving scales around quasars at z < 0.4 
(see also iHennawi etHI 120061: IMvers et al.] I2008D . For 
comparison, the median third-nearest-neighbor distance 
for our sample of early-type galaxies (over all environ- 
ments) is ~ 0.6 h^^ comoving Mpc in projection. In 
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addition, iCroom et al.l (|2005[ ) find that the average dark 
matter halo mass for quasars at intermediate redshift 
is roughly consistent with (within a factor of a few of) 
the minimum ha los mass inferred for groups a t z ^ 1 
(|Coil et alJ [20061 see also iHopkins et alJ l2007[ ). Alto- 
gether, clustering and environment studies do not sup- 
port (though also do not clearly exclude) quasars as a 
viable mechanism for size evolution at z < 2; regardless, 
questions still remain as to how quasar activity would 
be fueled in a massive early- type system at 2: > 1, as the 
standard scenario i nvolving the ma jor merger of two gas- 
rich s ystems fe.g.. iSpringel et alJl2005l : iDi Matteo et alJ 
[2005tlHopkins et alJ I2006ir fails to accurately describe 
massive early- type systems at z > 1 , which are relatively 
gas-poor and have stellar populatio ns with luminosity- 
weighted ages of > 1-2 Gyr (e g.. iDaddi et al.l 120051: 
Longhctti ct al. 2005; ' Treu et al.ll2005bt iSchiavon et al.l 
~006 : Combos ct al. .2007^ 

In contrast to quasars, Seyfert galaxies at z ~ 1 
typically reside in higher-density e nvirons, similar to 
those of gala x ies on the red seq u ence (IGeorgakakis et al.l 
20071 I2OOI iCoil et all [20091 iBradshaw c t al.l 120111: 
Digbv-North et al.l 12011 . but see also Silverman et all 



2009( 1. In addition, systems exhibiting line ratios consis- 



tent with Low Ionization N uclear Emission-line Regions 
fLINERS. [Heckman|[l980t) at z - 1, which tend to re- 
side on the red sequence, are likely to inhabit slightly 
more overdense regions even relative t o galaxies of like 
color and luminosity (i.e., stellar mass, lYan et al.1 120061 
[20Tlt iMontero-Dorta et al.l[2009t [juneau et al.ll201lD — 
though, recent work suggests t hat LINERs may no t 
be the product of AGN activity ()Yan fc Blantonll201lD . 
While lower-luminosity AGN are found in high-density 
regions at z 1, consistent with being the driv- 
ing mechanism behind the observed size evolution of 
early-type galaxies, the lack of variation in the out- 
flow velocities of winds observed in AGN hosts ver- 
sus star-forming galaxies suggests that low-luminosity 
AGN may not pla y a dominant role in galacti c feed- 
back (iRupke et al.ir20 05: Wcincr et al. 2009; Rub in et al 



201c 



20111 ; ICoil etal.ll2011L but see also iHainline et al 



20111 ). Furthermore, the mass loss needed to prod uce a 



factor of > 2 increase in siz e is of order 30%-50% ((Zhaol 
I2OO2I : iHopkins et al.|[20"l0aD , beyond the expected impact 
of fee dback from lower-lum inosity AGN or evolved stars 
(e.g.. lDamianov et al.l[2009[ ). 

5. SUMMARY 

Using data from the DEEP2 and DEEP3 Galaxy Red- 
shift Surveys, we have completed a detailed study of the 
relationship between galaxy structure and environment 
on the massive (10 < logio(M^,//i"2 Mq) < 11) end of 



the red sequence at intermediate redshift. Our princi- 
pal result is that at fixed stellar mass, redshift, Sersic 
index, and rest-frame color we find a significant rela- 
tionship between galaxy size and local galaxy density at 
z ^ 0.75, such that early- type galaxies in high-density re- 
gions are more extended than their counterparts in low- 
density environs. This result is robust to variations in 
the sample selection procedure, including selection lim- 
its based on axis ratio, surface brightness, and Sersic 
index. The observed correlation between size and envi- 
ronment is consistent with a scenario in which minor, 
dry mergers play a critical role in the structural evo- 
lution of massive, early-type galaxies at z < 2 and in 
which the evolution of massive ellipticals is accelerated 
in high-density regions. Future work, for example from 
the CANDELS HST/WFC3-IR imaging program, wiU 
soon enable complementary analyses at yet higher red- 
shift and in r nore extreme environm ents such as massive 
clusters (e.g.. iPapovich et"aI]|201lD . 
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